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The London penetration depth for the unconventional superconductor S^RuCU is analyzed 
assuming an order parameter which breaks time reversal symmetry and parity simultaneously. 
Such a superconducting state possesses chiral quasiparticle states with subgap energies at the sur- 
face. We show that these subgap states can give a significant contribution to the low-temperature 
behavior of the London penetration depth yielding a T 2 power-law even though bulk quasiparti- 
cle spectrum is gapped. The presence of several electron bands gives rise to interband transition 
among the subgap surface states and influences the properties of the surface impedance. Fur- 
thermore, the surface states lead also to a non-linear Meissner effect. 
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Several years of intense experimental research have es- 
tablished the unconventional nature of superconductivity 
in S^RuO/t. 1 ' 2 ) This compound has a layered perovskite 
structure representing a basically two-dimensional metal 
with three almost cylindrical Fermi surfaces. The sym- 
metry of the superconducting state is very likely odd in 
parity, which implies the spin-triplet configuration anal- 
ogous to superfluid 3 He. 3,4 ^ Muon spin rotation experi- 
ments provide evidence for broken time reversal symme- 
try, 5 ) a fact that strongly suggests that the gap function 
has the basic form, 



d(k) = A z 



(i) 



which is a chiral p-wave state, here written in the vec- 
tor representation, assuming cylindrical symmetry. The 
Cooper pairs possess an internal orbital angular momen- 
tum which is oriented along the z-axis. A consequence 
of this topological property of the superconducting phase 
is the presence of chiral surface states at the surface. 6 ~ 8 - ) 
While the chiral p-wave state has a basically gapful quasi- 
particle spectrum, the surface states correspond to sub- 
gap quasiparticle excitations with a continuous spectrum 
down to zero energy. 7, 8 ) These quasiparticle states are 
Andreev bound states and extend only over a coherence 
length towards the bulk. In this letter we consider the 
contribution of these states to the temperature depen- 
dence of the London penetration depth. The London 
penetration depth A|| for currents within the plane and 
the in-plane coherence length £m are very similar giving 
a Ginzburg-Landau parameter k — ~ 2.6. There- 

fore, the presence of the surface states can lead to a visi- 
ble reduction of the screening effect and could even dom- 
inate the low-temperature behavior An, in particular, if 
the bulk quasiparticle spectrum is gapped. We show here 
that power-law temperature dependence can result from 
the surface states, which is usually taken as an evidence 
for nodes in the bulk quasiparticle gap. 

The discussion of the London penetration depth re- 



quires a careful analysis of the current-current response 
to an external field which can be written in general as 



c 

4-7T 



Y.j dt 'j d\'K^(v,t-v',t')A v {v',t') 

(2) 

where only the transverse component of A^ir' , if) enters. 
The kernel K^(r,t;r',t') is obtained from the current- 
current correlation function. In our case this response 
consists of two contributions: the bulk part due to the 
continuum of quasiparticle states above the gap and the 
part due to the surface states. We consider from now on 
the specific case of a surface with normal vector along 
the x-axis and an external field parallel to the z-axis. 
Consequently we have to deal with the transverse vector 
potential and screening current along the y-axis. The 
relevant terms are then, 



]y{v,t) = -- 



dS 



2,5(4) ( 7 



n w (r;r')- t "'^ '° j A y (r') 

(3) 

where the integral runs of the four coordinates r' = 
(r', t') and Ao corresponds to the bare "London pene- 
tration depth" of the bulk regime which can be consid- 
ered as basically temperature-independent for very low 
temperatures. For this bulk part we take the local ap- 
proximation, while for the first part connected with the 
surface states nonlocality is important, as we will see 
below. 

The surface states can be easily described within 
the Bogolyubov-de Gennes formalism if we neglect self- 
consistency of the gap Ao which we choose to be con- 
stant everywhere inside the superconductor. For the sake 
of simplicity we assume also that the surface provides 
specular reflection of quasiparticles and the gap has no 
anisotropy on the Fermi surface. The electron band in 
our model has cylindrical symmetry and is represented 
by the parabolic form, = {(k^ + ky) — fc|}/2m neglect- 
ing any dispersion along the z-axis. Using these simplifi- 



2 



Takao MORINARI and Manfred SlGRIST 



cations the wave function of the subgap states localized 
at the surface is given by 



"k (r) 
Vk(r) 



Jk y y- 



«o sin(k x x) 



1 



(4) 

Since only states very close to Fermi surface are im- 
portant for the low-temperature properties the wave 
vector can be represented essentially as (k x ,k y ) = 
kF(cos6,sm0) for &f£o 1. Further, L y is the ex- 
tension of the system along the y-direction with peri- 
odic boundary conditions and d is the interlayer spacing 
(the wave function is renormalized per layer). The en- 
ergy of the surface states is given by E ky — r\A^k y jkF 
with r] = ±1 denoting the sign of the chirality of 
d(k) = z(k x ± ik y )/kF- 9 ^ Defining the current opera- 
tors as — (he/2mi)( 1 t> ji d ll ^ — h.c), where 



WW ~ I ^J(r) 



E 

k 



«k(r) 



<( r ) 



7kt 



(5) 

is the Nambu field operator with 7^ the Bogolyubov 
quasi-particle operator, we can express the current- 
current correlation function U yy (r]r') as 



n OT (r;r') 



TiV 
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lim 

>r, r'r' 



x(d Vl -d y2 ){d v , 2 -a Bi )Hr[G(ri;ri)G(/ 2 ;r 2 )] (6) 

where G(r;r') is the Nambu-Gor'kov Green's function 
in real space and d y denotes the derivative with respect 
to the spatial y-coordinate. The Green's function can be 
expressed as 

<j){x)(j>{x') ^ <7 



G(r,r';iw„) 



L y d 
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sE k 



( 7 ) 

with 0(a;) = ^/2/£ exp ^Z^ sinfci?a;, cto the unit ma- 
trix and &2 the second Pauli matrix, and u>„ being the 
fcrmionic Matsubara frequency. Using Eqs. (6) and 
(7), we calculate the current-current correlation function. 
The translational invariance along the y-direction allows 
us to transform the ^-coordinate into momentum space, 



ITyy (x, x , q^ zll n ) 



m 2 L y d 



g(x)g(x') 
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iCl n - E, 



k+q 



Eu 



^Ti 2 k% fk B T\ 2 g(x)g(x') 

in the limit T -C T c , where g(x) ~ exp(— |^) sin 2 (k F x)/£ 
is the square of the amplitude of the surface state wave 
function and iQ, n is the bosonic Matsubara frequency. 
In deriving Eq. (8), we restrict ourselves to the leading 
contribution for q <C kp and k B T <C Ao. The nonlocal 
nature of response enters via the product form g(x)g(x') 
which accounts for the fact that each of the quasiparticle 
state is localized at the surface. The field at the partic- 
ular point (x',y',z') couples to the surface state in the 
same layer with a weight g(x') and yields consequently 



a response at any other point (x, y, z') with weight g(x). 
This is a feature of the effectively one-dimensional char- 
acter of the surface states within each layer. A local 
approach in this place would underestimate the role of 
the surface states in the low-temperature response. 

Combining Eqs. (3) and (8), where we further use 
the analytic continuation iCl n — > fko + iS, with the 
Maxwell equation \7 2 A y (r,u) = — ^rjj,(r, a;) we obtain 
an integro-differential equation for A y (r,t). The bound- 
ary condition is given by d x A y (r,u)\ x = = B z (q,u>), 
where B z (q,ui) is the external magnetic field at the sur- 
face parallel to the z-axis. We solve this equation using 
an approximation g(x) ~ exp(— 2x/£o)/2£o in the inte- 
grand (we ignore the fast oscillations), which is certainly 
valid for &f£o 3> 1- This allows us to calculate the sur- 
face impedance, Z(q,w) = Ai:cE y {x = 0,q,u))/B(q,ui) = 
—4:7riu>A y (x = 0,q,u))/B(q,cj). We then obtain the pene- 
tration depth using the relation AttloX^q, lo) = lu\Z(q, u>). 
Taking a static limit Hukp/qAo — ► and q — > 0, we find 
in the regime of k B T <C A , AA(T) = A(T) - A(0) has 
T 2 -behavior: 

AMTyX^-^-^iksT/Aof. (9) 

Setting k ~ 2.6, which is a typical value of Sr 2 Ru04, 
we obtain AA(T)/A ~ 0.14 x (T/T c ) 2 , if we assume the 
weak-coupling relation A = 1.76fcsT c . We ignored the 
temperature dependence of Ao as it is exponential in the 
low-temperature regime in our model. 




Fig. 1. Schematic spectrum of the chiral surface states in the 
multi-band case of Sr2Ru04. The Brillouin zone contains three 
Fermi surfaces. The electron-like /3 and 7-Fermi surfaces yield 
a chiral surface state spectrum centered around k y = 0, while 
the spectrum of the surface state due to the hole-like a-Fermi 
surface located at the Brillouin zone boundary with opposite 
sign of chirality. The gap magnitudes are in general different. 
We assigned, however, the same magnitude for a and /3-Fermi 
surface for simplicity. 



It is important to notice that this contribution is in- 
dependent of the sign of the chirality and the charge of 
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the carriers, i.e., whether the superconducting state is derive the correlation function, 



d oc z(k x + iky) or z(k x — ik y )) and the Fermi surface 
is electron- or hole-like. Therefore, the formation of do- 
mains of the two superconducting states would not lead 
to a significant change of the result. Furthermore, for the 
case of several superconducting bands the contributions 
of the surface states of each band add up to enlarges 
the prefactor of the T 2 -law. The coherence length as 
the extension of the surface states towards the interior 
is different for each band, since the Fermi velocities and 
the gap magnitudes arc different. The coherence length, 
experimentally determined via the measurement of if C 2, 
giving k ~ 2.6 is the shortest among all. Therefore, the 
enhancement of the surface state contribution can be siz- 
able in the multiband case. We consider here the case 
of three bands, as schematically shown in Fig.l, where 
each band yields its own surface state described by a 
Green's function G^(r; r'), (the superscript j labels the 
j th ba,nd). We assume that the reflection of quasiparti- 
cles on the surface does not lead to transitions among the 
different bands. Each band is characterized by a Fermi 
vector kp\ the effective band mass rnj and the super- 
conducting gap Aj. We now use Eq.(6) to calculate the 
contribution of each band to the current-current corre- 
lation function and then analyze the resulting equation 
for the transverse vector potential as in the single band 
case. This leads to the low-temperature behavior of the 
London penetration depth, 



AA(T) 



3 ^ 



A A 



_ fk B TY 

3 ^ (2/sW) + l) 2 yX^Aj \ A J 



where 
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(10) 



(11) 



X /^> with ^> = Jl 2 k { p /mjAj. We 



and kW) = Xq/^P with 
choose A to reproduce again the weak coupling relation, 
Ao = 1.76fcsT c . The contribution of all three bands can 
easily give a prefactor to the (T/T c ) 2 -law of order one, 
consistent with recent measurements for fields along the 
z-axis. 

In the discussion of the multi-band situation we ne- 
glected the interband effects. Oscillatory fields, for ex- 
ample appearing in microwave experiments, can yield in- 
terband transitions. The matrix elements for the tran- 
sition depends on various details of the orbital and 
band structure. We do not go into these complex de- 
tails here, but assume that the interband transition can 
be described by an ordinary current operator, fy'^ — 
{he/2m'i)(i> < - i ^d^^ - h.c.) where m! is a phenomeno- 
logical parameter accounting for the matrix element. 

i s the Nambu field operator of the i-th band as in Eq. 
(5). In the small- momentum (q) and small- frequency (u) 
limit only surface states are important which share the 
same zero-energy momentum. As shown in Fig.l this is 
the case for the (3- and 7-band. The a-band is unimpor- 
tant because for interband transitions a large momentum 
transfer (q ~ ir) is necessary. Analogous to Eq.(8) we can 
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(12) 

exp(-2a;/|)/2 v /^ /3) ^ 7) with I" 1 = 
The surface state spectra in the two 



bands is approximated by E 



;k, with Vj = A 



We now take the analytic continuation iSl n — ► fiuj + iS 
and set q = in Eq.(12). For the limit of very small 
Huj (<C fcfiT) the surface impedance gets the following 
contribution from the interband transitions, 



ImZ ib (u>) = 



3<?C(7) 



4tt^ 



3A^ 2 (2/s + 1) 2 



A 7 m' 



2 - ( %L) In^l 

7T \A / V/3 



(13) 

for the imaginary part which yields in the zero frequency- 
limit also a T 2 -contribution to the London penetration 
depth (k = A/£ and v = v 7 — vp > 0). For the real part 
we obtain, 



ReZ ib (u)) 



2tt 3 |C (7) 



A^ 2 (2k +1) 2 



A m 7 \ ' 
A 7 m'7 



(14) 
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Juu\' 
A~o) 



Tiuo 
kWf 



for given T and small Two (-C fcgT -C Ao). The imag- 
inary part, the inductive resistance, shows a cj-linear 
plus w 3 -behavior, while the real part, the surface resis- 
tance follows an w 4 -law. In the opposite limit where 
ksT <C hu <C A the surface impedance due to inter- 
band transitions has to vanish. The reason is that for 
q = the initial and final states are either both empty 
or occupied in the zero-temperature limit. A simple anal- 
ysis shows the following behavior, 



ReZ ib (uj) oc — c~ 



2S ImZ ib (ui) oc -j (15) 



The surface resistance and the inductive resistance van- 
ish exponentially and with a power-law, respectively, in 
the zero-temperature limit. 

We now consider the possibility of a so-called nonlin- 
ear Meissner effect. The application of a magnetic field 
introduces a Doppler shift which changes the quasiparti- 
cle energy, 



EL = Et J) + —^-A v 



(16) 



if we again consider the case of n = (1, 0, 0). We can ex- 
pand the current-current correlation function for small 
A y and analyze the contribution to the London penetra- 
tion depth in the same way as done above. Restricting 
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to the single band model we obtain, 

AA(T) 4tt 2 k fk B T\ 2 ( 3k H \ 

(17) 

which leads to a non-linear correction in the external 
field. This effect is a consequence of the angular momen- 
tum of the Cooper pairs coupling to the field along the 
z-axis. The sign of this correction depends on the chiral- 
ity (direction of angular momentum) and the character 
of the Fermi surface, i.e., a different sign appears for the 
a-band than for the (3- and 7-band in Fig.l. Therefore, 
the presence of electron-like and hole-like Fermi surfaces 
as well as the formation of domains of the two chiral 
states lead to compensations which diminish the effect. 

In recent experiments by Bonalde et al., the temper- 
ature dependence of the London penetration depth was 
determined using a self-inductive technique. 10 ) It was 
found that the low-temperature behavior is indeed gov- 
erned by a T 2 -behavior. This is not compatible with a 
simple interpretation in terms of line nodes in the gap as 
originally proposed based on the T-power laws in specific 
heat and NQR, 11 ) since this would lead to a linear In- 
dependence. 12 "* The more sophisticated approach based 
on a nonlocal response theory by Kostin and Leggett 
(KL) (for the reason that k is small), however, would 
yield a T 2 -behavior. 13 - ) On the other hand, in this letter 
we propose an alternative mechanism for a T 2 -behavior 
based on the contributions of the surface states. In both 
theories it is expected that this power-law behavior is 
absent for A_l , for screening currents flowing along the z- 
axis. The z-axis current is not proportional to the surface 
state energy as required to obtain the T 2 -behavior in our 
theory. Furthermore, k± is about 20 times larger than 
the in-plane k so that contribution of the surface states 
as well as the nonlocal effect by KL are rather small. 
However, the measurements for fields in the plane, prob- 
ing the z-axis current show a similar T 2 -behavior. This 
is in apparent conflict with both interpretations. Since in 
this case, however, not only Aj_ but also the contribution 
from in-plane currents from the surfaces normal to the 
z-axis are involved, the final answer will be given only 
when these geometrical aspects have been thoroughly in- 
vestigated. 14 ^ 

In our model the gap size is isotropic on the Fermi 
surface. Anisotropy in turn modifies the surface state 
spectrum to Ek v — vk y + v'ky + ... without destroying 
the particle- hole symmetry (vv' < in general). This 
yields an additional T 4 -contribution which may not be 
so small. Together with other correction in this order 
this leads to 

^ffi = a(T/T c f + b(T/T c ) 4 + ... (18) 

with a ~ b > 0. In an intermediate range the second 
term generates a T-dependence which over some tem- 
perature range appears to be close to a T 3 -behavior and 
only at rather low temperatures the T 2 -law would dom- 
inate. For one sample Bonalde ct al. could indeed fit 
their data reasonably well with a T 3 -curve. 14 ' While 
this sample happened to be dirty it is not clear from 



the experiment what is the intrinsic origin for the appar- 
ently different power-law. Therefore, the additional T 4 - 
contribution, which due to the surface orientation and 
disorder is larger than usual, may be one possible expla- 
nation. 

A further aspect noteworthy here is that our mecha- 
nism is active at the surface only, while the KL scheme 
also applies also in the bulk of the superconductor. 
Therefore, the London penetration depth governing the 
magnetic interaction between vortices should have differ- 
ent temperature dependence in the two scenarios. Mea- 
surements of in the mixed phase by ^(SR suggest that 
London penetration depth saturates faster than T 2 at 
low temperature. 15 ) In contrast to the surface-sensitive 
experiment mentioned above, 10 -* ^(SR is indeed a bulk 
probe. Unfortunately, it has less accuracy in determin- 
ing the temperature dependence of the London penetra- 
tion depth so that we cannot draw a strong conclusion 
to date. Nevertheless, the present experimental result 
is consistent with the interpretation based on surface 
states. 
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